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This paper proposes a method to estimate state of charge (SoC) for Lithium- 
ion battery pack (LIB) with N series-connected cells. The cell’s model is 
represented by a second-order equivalent circuit model taking into account 
the measurement disturbances and the current sensor bias. By using two 
sigma point Kalman filters (SPKF), the SoC of cells in the pack is calculated 


by the sum of the pack’s average SoC estimated by the first SPKF and SoC 

differences estimated by the second SPKF. The advantage of this method is 
Keywords: the SoC estimation algorithm performed only two times instead of N times 
in each sampling time interval, so the computational burden is reduced. The 
test of the proposed SoC estimation algorithm for 7 samsung ICR18650 
Lithium-ion battery cells connected in series is implemented in the 
continuous charge and discharge scenario in one hour time. The estimated 
SoCs of the cells in the pack are quite accurate, the 3-sigma criterion of 
estimated SoC error distributions is 0.5%. 
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1. INTRODUCTION 

Lithium-ion batteries pack (LIB) are so far known to be the most commonly used batteries with the 
following advantages: a higher energy density, voltage capacity; lower self-discharge rate than other 
rechargeable batteries; the price is much cheaper than the LiPo battery. LIBs power almost every portable 
devices such as smartphone, laptop, and tablet sold today across the world. Now, LIBs are applied to not only 
electric cars, electric motorcycles, buses, trucks but also large-scale electrical energy storage. LIB is formed 
by many cells connected in series and parallel in order to get the desired voltage and power [1] as shown in 
Figure 1. The largest market for LIBs has traditionally been portable electronic devices, but in the near future 
there is an extensive growth in the demand for LIBs in transportation and electrical energy storages [2]. The 
state of charge (SoC) is an important input for balancing, energy and power calculations of LIB. SoC 
depends on the open-circuit-voltage (OCV) with the nonlinear relationship affected by the temperature, 
aging, and charge/discharge mode, [3]—[5]. The OCV depends on temperature and electrode particle surface 
concentration, but SoC depend on particle average concentration. In the practice, optimal safe operation 
range of SoC for cells is above 20% and below 90%. Figure 2 describes a typical relationship between SoC 
and OCV at -5 °C and 25 °C of temperature for a Lithium samsung cell. SoC is a parameter that is not 
measurable directly, so it is estimated by using the current, voltage, and temperature of the cell. SoC is 
needed to be estimated accurately in order to use for the SoC balance control problem. Controlling the 
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operation of the battery system in order to prolong the life of the battery cells, to operate safely is a topic that 
scientists are currently interested in applying LIBs in transportation and electrical energy storage. To estimate 
SoC for one cell, we need the model describing the dynamic behavior of the cell and the SoC estimation 
algorithm based on that model with real-time current, voltage, and temperature inputs. The accuracy of the 
estimated SoC depends on the model of the cell and the SoC estimation algorithms used [6]. 

The model of the cell: if the complex model is used (such as an electrochemical model), the 
estimated SoCs will be more accurate, but the computational burden and hardware costs will increase a lot 
when the number of cells in series N is large. If the equivalent-circuit model is used, the order of that model 
plays an important role. The computational complexity of the SoC estimation algorithm is simpler with the 
lower-order model but the cell dynamic behavior is not reflected accurately. There is higher accuracy in SoC 
estimation when the higher-order model is used, but the computational burden is increased [7]. The SoC 
estimation algorithm: Up to now, most of the SoC estimation algorithms are based on coulomb counting [8], 
Kalman filter (KF) [9], adaptive extended Kalman filter (EKF) [10], [11], and other [12], [13] improved 
algorithms. Several other algorithms such as the use of fuzzy models and neural networks are also applied 
[14], [15]. 


Figure 1. Inside the tesla model S battery pack 
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Figure 2. The typical relationship between SoC and OCV for a lithium samsung cell: (a) tested at -5 °C and 
(b) tested at 25 °C 


The SoC estimation for N cell connected in series is the challenge problem, there are some works 
related to this in literature [16], [17]. Normally, the SoC estimation algorithm is implemented repeatedly N 
times, if the number of cells N is large then this is really a computational burden for the battery management 
system [18]. Beside that, from the works [7], [12], [14], [17] and the actual testing we found that the 
interference of the measuring noises, current sensor bias are an obstacle to estimating SoC exactly. In 
practice, the current through the cell is usually measured by a shunt resistor with a voltage signal of mV, this 
voltage signal is more susceptible to interference than the voltage measuring signal at both ends of the cell 
(2.7 V to 4.2 V). The current sensor bias of the rectifier affected by temperature needs to be considered in the 
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current measurement. The KF theory assumes that all noises are zero mean, so an unknown current sensor 
bias can introduce permanent SoC error. This paper introduces the SoC estimation method for a LiB 
connected in series using sigma point Kalman filters (SPKF) to minimize the influence of current and 
vfoltage measurement noise on the estimated SoC results. By using the second-order equivalent circuit model 
(SECM), this allows to describe more exactly the cell’s dynamic in the charged and discharged scenarios 
with amplitude varying suddenly (an important feature of electrical vehicles using LIBs must be considered) 
[19]. This method uses two SPKFs, the SoC of cells in the pack is calculated by sum of pack’s average SoC 
estimated by first SPKF and SoC differences estimated by second SPKF. 

The contributions of this work are: The first, by using only two filters to estimate the N SoCs of N 
cells connected in series, the SoC estimation algorithm is performed only two times instead of N times in 
each sampling time interval, so the computational burden reduced to N — 2 times, this is a very effective 
solution when the battery pack has a large number of cells. The second, the SoC estimation taking into 
account the current sensor bias in the model in association with applying the SPKF filter allows mapping the 
distributions of current, voltage, and temperature signals to the distribution of the estimated SoCs of the cells 
with more accuracy tested by the 3o criterion. The remaining part of this paper is structured as follows: 
section 2 mentions the SECM of the cell taking into account the measurement noises and the current sensor 
bias. Section 3 introduces the proposed SoC estimation algorithm of cells connected in series using 02 
SPKFs. Section 4 shows the results applied on the Lithium ICR18650-22P SAMSUNG cells and discussions. 
Finally, section 5 summarizes some conclusions. 


2. THE SECOND-ORDER EQUIVALENT CIRCUIT MODEL OF CELL 

Suppose that the LIB’s cell is depicted as in Figure 3(a), in which i(t) is the cell current (i(t) > 0 
for charge and i(t) < 0 for discharge, i(t) = 0 for cell in no charge/discharge), 10% < z(t) < 95% is SoC, 
3.6V < v(t) < 4.2V is the cell voltage, —5°C < T(t) < 50°C is cell temperature. In order to estimate the 
SoC at the sample time k, which is z,, we need to measure the voltage, current and cell temperature. Let 
İk + by + nj, is the measured current, Vg + Ny is measured voltage, in which bẹ is the current sensor bias; 
Ni is the current noise; ny x is the voltage noise. The Figures 3(b), 3(c) are the recorded voltage and current 
in the practice, respectively. The scanning electrochemical microscopy (SECM) of the cell is plotted in 
Figure 4 with OCV(z(t),T) is the relationship between the oral cholera vaccine (OCV) and SoC, this 
relationship is nonlinear function depending on cell temperature and this function is needed to determine by 
the practice data collected by tests. The parameter h(t) is the dynamic hysteresis voltage as a function of SoC 
and time, R,(T),C,(T), R2(T), C,(T) are two resistor-capacitor pairs describe the cell dynamic, R(T) is 
series resistance parameter, Mj(T), M(T) are the hysteresis parameters, n(T) is Coulombic efficiency, y(T) 
is the parameter related to the the rate of decay. Note that all parameters of cell depend on the cell 
temperature. In exact terms, they also depend on the cell aging, the charge/discharge condition. For different 
kind of cell (Samsung, Panasonic), they are different. In this work, we suppose that all parameters of cell 
only depend on the cell temperature. 

Suppose that all the quantities is discreterized at the sampling time k,k = 0,1,2,...,00 with the 
sampling interval At is small enough so that the current and voltage of cell are constant over that interval. 
Suppose that Q is total capacity of a cell, the change in SoC over sampling interval is (1). 


. Ngåt 
Zkt = Zem k g (1) 


The diffused resistor currents are computed as (2), 


; = At), 1 At 3 
IR k+1 = CxXp\— RG. IRyk + — exp \— RG lk 


ire+1 = exp (— =) ir + (: — exp (- £) ix (2) 


Hysteresis voltage is given by (3). 


se) ) he = (: Š exp ( 


ices = exp (- sete) sign i) 3) 
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Figure 3. Illustration of the cell measurements, (a) the LIB, (b) the recorded voltage (c) current of cell 
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Figure 4. The second-order equivalent circuit model of the cell 


The current sensor bias is by, = bg-1 + Ng-1, in which Ng-1 is virtual noise of the current sensor 
bias used only for the calculation. The output equation of the model describing the relationship of the cell 
voltage, the current, SoC, the hysteresis, and the diffused resistor currents is written as (4). 


Yr = OCV (Zk) + Mhy — Riig, — Rir, — R2(ik — bg) (4) 


The state vector related to the diffused resistor currents, the hysteresis, SoC and the current sensor bias as is 
formed as the (5). 


Xp = [İr k+1 İr,k+ı hk+ Zk+ Desa]? (5) 


Define the system matrices and output matrix of the model as. 


E exp (--=) 0 E 1- exp (=) P 
Arc = i =e oi z = ‚Brc = on z , Ay, = exp ( aeut) 
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Suppose the input and output vectors of the model as (6). 


=| ik + bk 
k= 


sign(iz + al Yk = Vk (6) 


Finally, the SECM of the cell is written as (7). 


í Xk+1 = Alix) Xp + Blix) Up (7) 


Yk = OCV( Zk» Ty) + Cx, + Du, 


Model (7) is used to estimate the SoC for the cell with input ux and output Yg without measurement 
noise. When measuring Ug and Yg, the two measurement noises Nip and Ny, were mixed in the curent and 
voltage as described in Figure 3(a). So that, the input vector and output vector of the model is written as (8). 


İk + bk + Nik 
nk = > Yank = Ve tNy (8) 


sign(i, + by+Nix) 


With the input and output vectors in the (8), the model (7) is rewritten as (9). 


Xnk+1 = A(ink)Xnk + B(inx)Unk (9) 
Yn = OCV( Zn» Tk) + CXn e + Dunk 
In which x, +1 is state vector affected by measurement noise formed as (10). 

Xn = [İR n,k+1 lRonkt1 Ankst Znk+1 Bnk] (10) 


Suppose the error Wx = Xk+1 Xn,+1 to be the system noise of the model, § = Y,—Y,,, is the 
output noise of the model, we have (11). 
Xie1 = Alin) ¥Xnk + Bling) Une + We an 
Yk = OCV ( Zk, Tg) + Cx, + Dunx + čk 


The model (11) is used to estimate the state vector Xķ+1, in that vector there are the SoC Zę and the current 
sensor bias by need to be estimated. In this work, we use LIB Samsung ICR18650-22P [20] with parameters 
Q = 2200 mAh, Vnax = 4.25 V, Vmin = 2.75 V, nominal voltage 3.62 V. The test temperature range is 
considered from —5 °C to 45 °C. The values of model’s parameters determined by practical tests are given in 
Table 1. 


Table 1. The cell’s model parameters of LIB Samsung ICR18650-22P 


T°C -5°C 5°C 15°C 25°C 35°C 45 °C 
Nk 1.0869 0.9803 1.0220 0.95 1.0542 1.0399 

Q (Ah) 2.1596 2.1877 2.1943 2.175 2.1515 2.1523 

Y 243.1946 78.4915 63.6762 250.0074 170.6407 151.3064 

Mo (V) 0.0072 0.0049 0.0035 0.002 0.0036 0.0024 
M (V) 0.0347 0.0257 0.0218 0.02 0.0201 0.0185 
Ro (Q) 0.013 0.023 0.0292 0.033 0.0012 0.0011 
R,C, (S) 0.6124 1.7555 0.3227 0.899 0.6997 0.4630 
R,C, (S) 3.9035 7.5994 8.1118 10.023 8.1840 6.5319 
R, (Q) 0.0204 0.0203 0.0201 0.019 0.019 0.0019 
R2 (Q) 0.0494 0.0376 0.0288 0.022 0.0136 0.0134 
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3. PROPOSED SOC ESTIMATION ALGORITHM FOR N CELLS CONNECTED IN SERIES 
USING TWO SPKFS 

Suppose that there is N cells LIB connected in series as shown in Figure 5, the problem is to 
estimate SoCs of N cells with the input and output signals affected by noises as (8), the noise distributions 
are considered as Gaussian-distributed. SPKF is the filter with high precision aimed at nonlinear system, in 
practice, the SPKF has some advantages compared to the EKF. This filter allows performance KF for 
nonlinear systems without the linearization steps. The set of weighted samples, called sigma points, is used 
for normalizing mean and covariance of a probability distribution. The sigma points of input distributions are 
mapped thought the nonlinear dynamic model functions to get the sigma points of output distribution 
[21], [22]. SPKF allows estimation states of nonlinear dynamic systems more accurately than the EKF does, 
but the computational complexity is same as EKF [23]—[25]. SPKF is easily implemented because there is no 
requirement of the Jacobi matrix calculation. The accuracy of state estimation using the SPKF mainly 
depends on the strategies of choosing sigma-points [26], [27]. Normally, we need to estimated SoC N times 
for N cells. If the number of cells N is large, there is a large amount of calculation [18], [19]. To reduce the 
SoC estimation times, in this work we use only two SPKFs, the first filter estimates the mean of N SoC and 
the current sensor bias, the second filter estimates N SoC differences from the average SoC, the estimated 
SoC of each cell then is calculated by the sum of the average SoC determined by the first filter and the 
difference SoC estimated by the second filter. In one sampling time interval, the SoC estimation algorithm is 
implemented two times instead of N times. 


1 
1 
1 Uy 
| ————> 
I 
Cell Cell Cell Cell 
al, 2 3 N 


Figure 5. N LIB cells connected in series 


Consider the second model of cell in the (11) with input and output in the (8). The SOC estimation 
algorithm flowchart for N cell connected in series is conducted as: 


The SOC estimation algorithm flowchart 
Begin 
Initialization 
For each sampling time k,k =1,2,3,...,00, do 
Measure the current of cell ing ; Measure voltages of N cells Yar Vets Vie ; 
Measure ambient temperature of cell Tk 
Estimate state vector Re = [Êr e+ Îr, k+1 Rear 2ra Benl by first SPKF 
For each cell i=1,2,3,..,N do 
Estimate the SOC difference AZO by second SPKF; Estimate 
the SoC for cell 20 =4, + a20 
End 


End 


3.1. Algorithm to estimate the average SoC for the N cells by the first SPKF filter 
The first filter estimates the mean of N SoC and the current sensor bias as steps following: 
— Form the augmented state vector at sampling time k defined by Rt, including the estimated state 


vector at the sampling time k — 1 the mean vector of the system noise W and mesurement noise ¥: 


get = Ria w E (12) 


— Form the augmented covariance matrix (CM) of state estimation error vector by using CM of state 
estimation error at sampling time k — 1, the CM of system noises and the CM of measurement noises: 
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Ogk-1 = diag|OZe-1 Sw 9z] (13) 
— Set up matrix of p + 1 sigma points: 
A A A A T 
xe = faat gat +y foot ety lott a Rn wA xt] (14) 
— Calculate sigma points matrix at sampling time k by using the cell model’s dynamic function in (11): 
Xe = Angling) Xe + Belin )UnktXen; =f = 0,1,2,..p (15) 
— Estimate the priori state vector at the sampling time k: 
Re = Dyno XG (16) 


— Update the CM of the state estimation errors: 


= Z VAS Zo Rone 
Oze = Djo % (Xx — Re ) KKJ — Ke ) (17) 
— Calculate the output sigma points matrix at the sampling time k by using the output equation in the 
model (11): 
Ye = OCV (Bp, Ty) + ChXT + Dating, + XE, (18) 


— Estimate the output voltage of cells as (19): 
Pk = Dio j” Yk (19) 

— Update the CM of output voltage estimation errors: 
OFK = Xo af (Yk — Ir Yk — Pr)” (20) 

— Update the covariance matrices of the state vector estimation error and output voltage estimation error: 
Ox ge = Lino % (Xj — Re) Me = Fe)” (21) 

— Update the estimation gain matrix as (22): 


Ly, = Z2 (22) 


Fk 


— Update the estimation state vector by using the estimation gain matrix, the measured real output 
voltages and the estimated output voltage of cells: 


Re = Ke tLe Ve - Ik) (23) 
— Update the CM of state estimation error vector: 
Okr = ORK — Le OF scbk (24) 


3.2. SoC differences estimation for N cells by second SPKF 
The second SPKF estimates the SoC differences AZ, i=1,2,...,N, for N cells at the sampling 


time k as following steps: 


Int J Elec & Comp Eng, Vol. 12, No. 2, April 2022: 1334-1349 


Int J Elec & Comp Eng ISSN: 2088-8708 O 1341 


go 


— Form the augmented state vector of cell i at the sampling time k — 1 by using the SoC difference AZ 


and the mean vector of the SoC difference noises wÊ.: ), 


ADO = AEP w (25) 


— From the augmented CM of SoC difference estimation error of cell i, by using the CM of SoC 
difference estimation error of cell i at sampling time k — 1 and the CM of SoC difference noise: 


a,(i)+ 


St Oe 
T = diag (ea 0,0 ) (26) 


— Set up the matrix of p + 1 sigma points of cell i: 


T 
^ A sa, (i (i A Az,(i 
Axe = (AR Age O* +y Jose Aaa Y oka £ (Axes Ane i} eD 


— Calculate the state sigma point matrix of cell i at the sampling time k: 


A -i axe 
artO- = an ic addi ae 


AX j =0,1,2,..,p (28) 


3600Q 


— Estimate the SoC difference of cell i at the sampling time: 


Ag =7f o al Ay ce (29) 


— Update the CM of SoC difference estimation error of cell: 
S apa > GBAT 
Chea = yA oai af (Axi? i)— — az ) (An? = AZ ) (30) 
— Calculate the output sigma points matrix of cell i: 


YO = OCV (Zp + 420, Te) — inp ARË + 
MOR, — RË CM ig, — RP CP ing — Roling — tf) (31) 


— Estimate the output voltage of cell i: 


pË = ae ye (32) 


— Update the CM of output voltage estimation error of cell i: 


i ; . SAT 
aye = Loaf (Ve -IPP - 2P) (33) 
— Update the CM of SoC difference estimation error and output voltage estimation error of cell i: 
= D- Do LOT 

alya = Lian Hf A — A?” ) (Ye? — 9°) (34) 

— Update the state estimation matrix of cell: 
@ 
(i) _ [Az,5,k 
Le = (35) 


— Estimate the SoC difference of: 
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Mat = a2 + LO VK) (36) 
— Update the CM of SoC difference estimation error of cell i: 
; : \ Gy) .at 
Ose = Oker — Lye Oyy Lye (37) 


3.3. SoC estimation of cells 
The SoC of cell i calculated as (38): 


20t = at + ng* (38) 


4. RESULTS AND DISCUSSION 

The experimental system of SoC estimation is described in Figure 6, the components in the 
experimental system are listed as: i) The LiB samsung is formed by 7 cells connected in series with 
parameters listed in Table 1; measurement modul is created by 7 voltage difference measurement circuits, 1 
current measurement curcuit, 1 temperature sensor. The voltage difference signals are to the digital signals 
10 bit and transferred to the personal computer (PC) by RS485; ii) the PC implements the SoC estimation 
program for cells, the configuration of PC is core (TM) i3-6100 CPU@3.70 GHZ, 4 cores, RAM 8 GB; the 
MATLAB software is used to implement the SPKF algorithms, the LabVIEW 2018 [22] is used to make the 
graphic user interface, the MATLAB software is embedded to the LabVIEW to run the SoC estimation; 
iii) the load is motor BLDC, YONG-d67 240 W, 24 V, 10 A, the maximum speed is 295 rmp; iv) the charge 
source is owon, 195 W with 2 programmable channel 36 V, maximum current 6 A; and v) the SoC estimation 
algorithms are coded in MATLAB script, then embedded to the LabVIEW environment. The sampling time 
interval is 1 second. The algorithm flowchart of SoC estimation for N cell in shown in Figure 7. 


Implemented in the 
measurement modul and 
LabVIEW software 
Implemented in the Matlab software Implemented in 01 
fe cuieenuwodedcuuceseuewencecacewunewosuuar sampling time 
Implemented in the Lab VIEW software interval 


Figure 7. The algorithm flowchart of SoC estimation 


Int J Elec & Comp Eng, Vol. 12, No. 2, April 2022: 1334-1349 


Int J Elec & Comp Eng ISSN: 2088-8708 O 1343 


4.1. The SoC estimation results for 1 cell 

In this part, the SoC estimation results for 1 cell is presented, the Figure 8(a) describes the current, 
voltage and temperature of the cell. The noise distributions of measured currrent, voltage and temperature are 
plotted in Figure 8(b), the Figure 8(c) and the Figure 8(d). The current noise was distributed inside from 
—0.8 A to 0.8 A caused by shunt-resistor with the sensitivity 15 mV/A. The range of the voltage noise is 
range of -0.03 V to 0.03 V, the distribution of temperature noise is in the range from —0.7 °C to 0.7 °C. 
Figure 9 shows the results obtained when performing the SoC estimation algorithm corresponding to current, 
temperature and voltage inputs as shown in Figure 8 using the second-order model of cell with the 
parameters described in Table 1. 

The Figure 9(a) describes the comparison between real cell voltage and estimated voltage calculated 
by cell model, the voltage error is plotted in Figure 9(b) with 0.1 V and —0.075 V of maximum for discharge 
and charge, respectively. This maximum error is caused when there is the changing of charge to discharge 
vice versa. The cause of voltage error is the dependence of model on the temperature varying in the test, the 
parameters of cell were determined at particular temperatures —5°C,0 °C,5 °C,10°C,15 °C, 20°C, 
25 °C, 30 °C, 35°C. In the test, the temperature of cell raises from 25.5°C at to 31.5°C as plotted in 
Figure 9(c). The Figure 9(c) also plotted the estimated SoC of cell. After 1 hour with the discharging (3 A) 
and the charging (1.5 A) alternately, The SoC of cell tends to decrease from 92.5% at beginning to 65% at 
the end of the test. The SoC estimation error is shown in Figure 9(d) by comparing to criteria 2, + 30. 

The criteria 30 decreases gradually from the beginning of the test to the end of the test as plotted in 
Figure 9(e), we can see that the accuracy of SoC estimation of the cell described by the criteria 30 at the end 
of test was 0.0125%. The Figure 9(f) plots the updating of matrix ož by the time. At the times, there are the 
changings of charge to discharge vice versa, the voltage errors are large, some elements of matrix oj need to 
be changed in order to keep estimating SoC smoothly. 


Voltage (V) 50000 
Voltage B Temperature 
574 AW) 437 CON 
4 \ Current (A) d H32 
N 40000 
H 
2 
Q 
© 4 30000 
2 z 5 
= 5 S 
g $ o] 
2 5 2 20000 
E R E 
z £ H] 
8] 5 2 
= 
10000 
t T T Too —r r T T T 0 
0 500 1000 1500 2000 2500 3000 3500 - 
Time (s) Distribution of current noise (A) 
(a) (b) 
50000 
25000 + 
40000 + 
20000 4 
g g 
g 3 
E 300004 5 
5 % 150004 
5 ja 
g 3 
2 4 
2 20000 — 10000 
3 
3 Žž 
100004 50004 
04 n 0 i r 
1 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 
Distribution of temprature noise (°C) Distribution of voltage noise (V) 
(c) (d) 


Figure 8. The measurements and distributions, (a) the current, voltage and temperature of cell in the test, 
(b) the noise distributions of current, (c) temperature, and (d) voltage 
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Figure 9. The SoC estimation results for 1 cell, (a) the voltage and estimated voltage, (b) voltage error, (c) the 
estimated SoC and temperature, (d) estimated SoC+36, (e) criteria 30, , and (f) the matrix Ofk 
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4.2. The SoC estimation results for 7 cells 

In this part, the SoC estimation results for 7 cells connected in series are presented. Suppose that all 
cells are the same with parameters listed in Table 1. The Figure 10(a) plots the cells current, the Figure 10(b) 
and Figure 10(c) depicts the voltages of cells and temperature, respectively, the distributions of noises in 
current, voltages and temperature are plotted in Figures 10(d), 10(e) and 10(f). The estimated SoC for 7 cells 
is presented in Figure 11 and Figure 12 (in appendix), in which the Figure 11(a) is the estimated average SoC 
for 7 cells, the Figure 11(b) plots the estimated SoC differences of cells. 
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Figure 10. The measurements and distributions, (a) the current of 7 cells, (b) the voltages of cells, (c) the 
temperature of cell, (d) the distributions of noise in current, (e) voltages and (f) temperature are plotted 


The Figure 12 plots the estimated SoC of each cell and comparing with criteria 30, in which: 
Figure 12(a) and Figure 12(b) are the estimated SoC and comparing with criteria 30 for cell 1, similarity for 
the cells from 2 to 7 are plotted in Figure 12(c) to Figure 12(0). The Figure 13(a) shows the varying of matrix 
až in the test, the Figure 13(b) plots the varying of criteria 30. The mean of the criteria 30 for cells at the 
end of the test is approximately 0.5%. Figures 12 and 13 see in appendix. 
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Figure 11. The estimated SoC for 7 cells, (a) the estimated average SoC, and (b) the estimated SoC difference 
of cells 


5. CONCLUSION 

This paper deals with the method for SoC estimation of LIB with cells connected in series. In this 
works, the cell’s model is represented by the second-order equivalent circuit model taking into account the 
measurement disturbances and the current sensor bias. Using 02 SPKF filters to estimate the average SoC of 
cells, the current sensor bias and the SoC differences, the estimated SoC of each cell is is calculated by sum 
of pack’s everage SoC estimated by first SPKF and SoC differences estimated by second SPKF. This allows 
reducing the SoC estimation times to 2 times, this is a very effective solution when the battery pack has a 
large number of cells. The practical test is applied to 7 samsung ICR18650-22P lithium-Ion battery cells with 
the discharging (3 A) and the charging (1.5 A) alternately in 1 hour. The estimated SoC is quite accuracy 
with the criteria 3 o of estimated SoC error is approximately 0.5. In this work, we suppose that all cells are 
physically identical, the equivalent circuit model with parameters listed in Table 1 are linearized at some 
particular temperatures and applied for all cells. This is clear evidence that those factors contribute to the SoC 
and current sensor bias estimation error observed in the test. In future work, the model of each cell to be 
updated adaptively by temperature will be researched. 
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Figure 12. The estimated SoC for each cell and the comparing with criteria 30, (a) the measurements for cell 
1, (b) estimated SoC for cell 1, (c) the measurements for cell 2, (d) estimated SoC for cell 2, (e) the 
measurements for cell 3, and (f) estimated SoC for cell 3, (g) the measurements for cell 4, (h) estimated SoC 
for cell 4, (i) the measurements for cell 5, (j) estimated SoC for cell 5, (k) the measurements for cell 6, (1) 
estimated SoC for cell 6, (m) the measurements for cell 7, and (n) estimated SoC for cell 7 
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